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ABSTRACT 

We present a unified model of optical and ultraviolet (UV) light curves for one of the fastest classical 
novae, V838 Herculis (Nova Herculis 1991), and estimate its white dwarf (WD) mass. Based on an 
optically thick wind theory of nova outbursts, we model the optical light curves with free-free emission 
and the UV 1455 A light curves with blackbody emission. Our models of 1.35 ± 0.02 Mq WD 
reproduce simultaneously the optical and UV 1455 A observations. The mass lost by the wind is 
A A/wind 2 X 10^^ Mq. We provide new determinations of the reddening, E{B — V) = 0.53 ± 0.05, 
and of the distance, 2.7 ± 0.5 kpc. 

Subject headings: — nova, cataclysmic variables — stars: individual (V838 Herculis) — stars: mass 
loss — ultraviolet: stars — white dwarfs 
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INTRODUCTION 
V83 8 Her was discovered independently by ISugano I 



from the Palom ar Sky Surve y plates as 19 and 17.5, 
respectively, by iWest'l (I1991D . and 20.6 and 18.25 by 



(119911) o n 199 1 March 24.78 UT at 5.4 mag and by 
lAlcock I ([1991 at ~ 5 mag on 1991 March 25.67 UT. 
Alcock's visual discovery around the peak was made in 
strong twilight, so the peak magnitude was not accu- 
rate. The outburst time can be estimated fro m the 
upper limit (9 mag) prediscovery observation by lUeta I 
(1991') on JD 2,448,339.9 and Sugano's discovery on JD 
2,448,340.281. In the present work we adopt the out- 
burst time JD 2,448,340.0 as day zero and the maximum 
magnitude my, max = 5.4 mag. 

Soon after the optical maximum, the nova entered a 
very rapid decline phase, followed by a slight oscilla- 
tory behavior enduring several days and, later on, by 
a smooth decline. Infrared (IR) fluxes also showed a 
rapid decline, reaching a local minimum on day 6, fol- 
lowed by a rapid brightening. The spectral energy dis- 
tribution was consistent with free-free emission from 
day 1.3 to 6.5, and thereafter with blackbody emis- 
sion, ascribed to the forma ti on of a hot dust shell 
([Chandrasek har et al.l [19921: iHarrison fc Stringfellow 
1994; Kidger fc Martinez-Roger II 1993t I Woodward et al. 
19921) . Indication of IR emission from silicate 



was rep orted only a t a later phas e (Lvn ch et al. 



grams 
1 ll99l 



ISmith et al.'lll995^ . ISmith et al.~l (|1995l) suggested that 
silicate emission was due to a light echo by cold silicate 
grains deposited in a previous nova eruption. 
Red and blue preoutburst magnitudes were estimated 
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iHumphrevs et al.l (|199lD . but were contaminated by a 
closeby (~ 1") star of similar magnitude. V838 Her re- 
turned to the preoutburst ma gnitude V = 19.0 on da y 
403, and V = 19.2 on day 572 ijSzkodv fc Ingram l [l99?l. 

In this series of papers, we have applied our univer- 
sal decline law to many n ovae, including moderate ly fast 
novae (e.g., V1668 Cyg in lHachisu fc Kato I (|2006f ) [here- 
after referred a s Pape r I]) and a slow nova (GQ Mus in 
iHachisu et al. I (2008) [ Paper HI]), and have succeeded 
in reproducing simultaneously optical, infrared, ultravio- 
let (UV) 1455 A, and supersoft X-ray light curves. V838 
Her is, however, one of the fastest classical novae, so it 
is interesting to see if the universal decline law also ap- 
plied to such an extreme object. The main observational 
features of V838 Her are summarized in Tables [T] and |2l 
Section |2| reviews our observational results based on the 
lUE spectra. Our model light curves are briefly intro- 
duced in Section |3l The results of light curve fitting are 
given in Sectional Discussions and conclusions follow in 
Sections [5] and |6l In Appendix we give several examples 
of nova light curves, which show high degree of homology, 
as predicted by our models. 

2. UV OBSERVATIONS 

V838 Her has been monitored by lUE from day 2 to 
day 831, mainly at low resolution. A g a llery o f UV 



spectra can be found i n ICassatella et al. I (|2004a[ ) and 



IVanlandinghai n~et al. I (119961 ) 

In the following we revisit the problem of color excess 
E{B — V) of V838 Her, and describe the long term evo- 
lution of UV continuum and of emission lines. The UV 
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TABLE 1 

Observational Properties of V838 Herculis 



subject data reference 



discovery 


... JD 2,448,340.28 


Sueano (1991) 


nova speed class 


very fast 


IR minimum 


... JD 2,448,345.94 (day 6) 




t3 


5.3-6 days 


Harrison & Stringfellow (1994) 




5.4 mag 


Sugano fl991) 


A4'v,max from ta 


... -9.94 ~ -9.80 mag 


equation JSj 


distance 


... 2.8-10 kpc 


Vanlandineham et al. fl996): Harrison & Strinefellow Q99^ 


distance 


... 2.7 ±0.5 kpc 


this work 


FWHM of UV 1455 A 


4.7 days 


this work 


E{B-V) 


... 0.3-0.6 


Vanlandingham et al. fl996): Harrison & Stringfellow (199^ 


E(B~V) 


... 0.53 ±0.05 


this work 


dust 


very thin 


see ^Ol 


orbital period 


... 7.14 hr 


Ineram et al. ("1992"); Leibowitz et al. (19921 
Szkodv & Hoarcl7T994) 


H-burning phase 


... < 1 yr 



TABLE 2 
Chemical Abundance by Weight 



object X Y CNO Ne Na - Fe reference 

Sun 0.7068 0.274 0.014 0.0018 0.0034 Grevesse & Anders (1989) 

V838 Her 0.78 0.10 0.041 0.081 0.003 Vanlandingham et al. (1996) 

V838 Her 0.59 0.31 0.030 0.067 0.003 Vanlandingham et al. (1997) 

V838 Her 0.562 0.314 0.038 0.070 0.015 Schwarz et al. (2007 )' 



spectra were retrieved from the lUE archive through the 
INES {lUE Newly Extracted Spectra) system-^, which 
also provides full details of the observations. The use 
of the lUE INES data is particularly important for the 
determination of reddening correction because of the im- 
plementation of upgraded spectral extraction and flux 
calibration procedures compared to the previously pub- 
lished UV spectra. 

2.1. Reddening Correction 

As summarized by IVanlandingham et al. I ()1996D , the 
color excess of V838 Her has been determined by several 
authors usi ng different methods , based on the Balmer 
decre ment (jlngram et al. I Il992t IVanlandingham et al. I 
1 1 9961). the equ ivalent width of the Na I interstellar lines 
( Lvnch et al. lfT992i) . the ratio of the UV flux above and 
below 2000 A ()Starrfield et al. ||1992D and the assumed 
intrinsic color at maximum light ([Woodward et aL~l 
[199^ . The mean value of these determinations, assum- 
ing to have the same weight, is E{B — V)= 0.50 ± 0.12. 
Such a large error would cause large uncertainties on the 
distance and hence on the intrinsic parameters of the 
nova, as derived through a comparison with our models. 
This has driven us to revisit the problem of the reddening 
using an independent method that we consider more re- 
liable than the previous ones, because of relying directly 
on its best detectable effect, i.e. the strength and shape 
of the 2175 A broad dust absorpt ion band . 

The Galactic extinction curve (ISeaton I Il979( ) shows 
a pronounced broad maximum around 2175 A due to 
dust absorption. Since it takes the same value ^(A) = 
A{X)/E{B - V^) « 8 at A 1512, 1878, and 2386 A, the 
slope of the straight line passing through the continuum 
points at these wavelengths is insensitive to E{B — V) 
in a (A, log Fa) plot. This circumstance can be used to 

^ http://sdc.laeff.inta.es/ines/ 



get a reliable estimate of E{B — V) as that in which the 
stellar continuum becomes closely linear in the 1512-2386 
A region, and passes through the continuum points at 
the above wavelengths. From 8 pairs of short and long 
wavelength lUE spectra taken from day 6 to day 69, i.e. 
during or close to the nebular phase, we have in this way 
found E{B -V) = 0.53 ± 0.05, a value that we adopt 
in the following. Examples of lUE spectra of V838 Her 
corrected with E{B — V) = 0.53 are reported in Figure 
[T]for day 10, 12, 15 and 22. 

In principle, we could have estimated the E{B — V) 
color excess also from the observed relative intensities of 
the He II 1640 A Balmer Hue, and the 2734 and 3203 
A Paschen recombination lines, compared with theoret- 
ical ratios, as done, for example in the case of GQ Mus 
(Paper III). Such lines were however too faint to be mea- 
sured accurately. 

2.2. Evolution of the UV Continuum 

We have measured the mean flux in two narrow bands 
20 A wide centered at 1455 A and 2855 A, selected 
to best represent the U V continuum because little af- 
fected by emission lines (jCassatella et al. II2002D . Figure 
Oshows the time evolution of the F(1455 A) and F(2885 
A) fluxes and of the UV color index C(1455 - 2885) = 
-2.51og[F(1455 A)/F(2885 A)]. The measurements were 
made on well exposed low resolution large aperture spec- 
tra. Figure reports also, for comparison, the visual 
light curve obtained from the Fine Error Sensor (FES) 
counts on board lUE, once corr ected for the time depen- 
dent sensitivity degradation (see lCassatella et al. Il2004al . 
for details on the FES calibration). 

It appears from Figure[5]that flux maximum is reached 
later in the UV 1455 A band than in the 2885 A and 
visual bands, as a result of the progressive shift of the 
maximum emissivity towards shorter wavelengths. The 
progressive hardening of the UV spectrum stops after 



Light Curve Analysis of V838 Her 



3 



V838 Her E(B-V) = 0^53 




1 500 



2000 2500 
Wavelength (A) 



3000 



Fig. 1. — lUE spectra of V838 Her obtained at different dates 
(days after the outburst). The spectra have been corrected for red- 
dening using E{B — V) = 0.53. The vertical dotted lines represent 
the wavelengths AA 1512, 1878, and 2386 A at which the extinction 
law takes the same value. With the adopted value of reddening, 
the stellar continuum underlying the many emission lines is well 
represented by a straight line all over the full spectral range. Sat- 
urated data points in the emission lines are labeled with pluses. 

day 10, when the C(1455 - 2885) color index stabilizes 
around a value of zero. Both these fe atures are a com- 
mon p roperty of the novae studied by ICassatella et al. I 
((200l . 

2.3. Evolution of the UV emission Lines 

Figure [3] reports, as a function of time, the observed 
flux in the most prominent UV emission lines, as mea- 
sured by us from the available lUE low resolution spec- 
tra. The figure shows that the maximum emission pro- 
gressively shifts from low to high ionization lines, such 
that Mg II 2800 A is the first to reach maximum, fol- 
lowed by C II 1335 A, O I 1300 A, C III] 1909 A, N IV] 
1487 A, and N V 1240 A maximum. The time delay of 
maximum emission At with respect to the visual maxi- 
mum does in fact increase with increasing the line ioniza- 
tion potential x [eV] , as better shown in Figure |4] where 
At, once normalized to the ^3 time, is plotted as a func- 
tion of X for the different emission lines. Here, ^3 is the 
time in which the visual magnitude drops by 3 magni- 
tudes from the optical peak. Such a behavior is consis- 
tent with that observe d in a sample of s even CO novae, 
whose mean data from ICassatella et al. I (12005) are over- 
plotted for comparison. In particular, note that the max- 
imum emission in C III] 1909 A and N III] 1750 A takes 
place after a time of At/t^ « 2, as in CO novae. This 
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Fig. 2. — Evolution of the continuum fluxes at 1455 A and 2885 
A, of the ultraviolet color index C(1455 — 2885), and of the VpES 
visual flux of V838 Her. Fluxes are in units of 10~^^ erg cm~^ s~^ 
A~^, not corrected for reddening. Only the color index has been 
corrected for reddening using E{B — V) = 0.53. 

time has been identified in Cassatella et al. as the start 
of the pre-nebular phase. 

2.4. High Resolution Emission Line Profiles 

Good quality LUE high resolution spectra of V838 Her 
are available only at two dates: day 2.5 (LWP 19993) 
and day 15 (SWP41317). The most prominent feature 
in the long wavelength spectrum of day 2.5 is the broad 
P Cygni profile of the Mg II 2800 A doublet (Figure 
[5]) . A comparison with theoretical P Cygni profiles com- 
puted with the SEI (Soblolcv Exact Integration) met hod 
(|Lamers et al. Iil987: Groenewegen fc Lamers Ill989[ ) in- 
dicates that the terminal velocity of the wind is ~ 3000 
km s~^. The huge emission compared with the absorp- 
tion component, usual in novae, is due to efficient pop- 
ulation of the upper levels of the doublet by electron 
collisions. 

Another interesting line profile is that of the C III] 
1906.68-1908.73 A doublet shown in Figured which is 
the strongest emission line observed in the short wave- 
length spectrum of day 15. As it appears from the fig- 
ure, the C III] emission has a very complex profile. One 
way to interpret this profile is that the emitting region 
is very inhomogeneous in terms of density and velocity 
structure. This possibility cannot be ruled out a priori 
also in view of the presence of complex emi ssion line pro- 
files, e.g., in the hydrogen Balmer lines fsee lliiima et al. I 
l2009f ). The other possibility is that the C III] emission 
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V 838 Her 1991 = 5 days) 



V838 Her 




3340 8350 8360 8370 
JD - 2440000 



8380 8390 



Fig. 3. — Evolution of the observed fluxes in the most prominent 
emission hnes of V838 Her, in units of 10~^^ erg cm~^ s~^. Open 
dots in the N IV] and N V panels refer to features not due to these 
transitions. 
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Fig. 4. — Delay of the flux maximum of the emission lines, 
normalized to the time, as a function of the corresponding ion- 
ization potential for V838 Her (filled circles). The emission lines 
considered, in order of increasing ionization potential, are: Mg II 
2800 A, C II 1335 A, O I 1300 A, C III] 1909 A, N III] 1750 A, 
N IV] 1487 A, and N V 1240 A. We also show, for compari s on, th e 
average values of CO novae (crosses) in ICassatella et al. I ()2005I V 
In this case, the x-axis has been shifted by 0.5 eV to make the error 
bars visible and the O III 1660 A line is added between the N III] 
and N IV] lines. 



- LWP19993 
4 1 (day 2.5 




2760 2780 2800 2820 2840 
Wavelength (A) 

Fig. 5.— P Cygni profile of the Mg II 2800 A doublet from 
the lUE spectrum LWP19993 of V838 Her taken on day 2.5. 
The observed profiles have been fitted with a theoretical P Cygni 
profile, obtained through the SEI method jLamers et al. I [19871 : 
IGroenewegen fc Lamers1ll989l ). indicated as a thick line. The de- 
rived wind terminal velocity is ^ 3000 km s~^. Fluxes are nor- 
malized to the local continuum. Note the presence of the nar- 
row interstellar components of the doublet at about the laboratory 
wavelength. 

is eroded by overlying absorption from an outer shell. 
A likely source of overlying absorption is the principal 
and diffuse absorption components from the Fe III UV 
34 triplet at radial velocity of —500 km s^^ and —1600 
km s~^, respectively. Such components have already 
been confidently id entified to affect the C III] emission 
line in V1974 Cyg (|Cassatella et al. Il2004bn . 

3. THE MODEL OF NOVA LIGHT CURVES 

In this and previous papers, we have presented a uni- 
fied model for the IR, optical, UV, and supersoft X-ray 
light curves of several classical novae. Our models are 
based on the optically thick wind theor y of n ova out- 
bursts as described in iKato fc Hachisu I (|1994l ) and in 
Paper I. 

3.1. Optically Thick Wind Model 

After a thermonuclear runaway sets in, the photo- 
sphere of the white dwarf (WD) envelope expands greatly 
to a giant size with i?ph > 100 Rq. The envelope set- 
tles into steady-state around the optical peak. The fol- 
lowing decay phase can be described by a sequen ce of 
steady state solutions (e.g.. lKato fc Hachisu 11199^ . Us- 
ing the same method and numerical techniques as in 
IKato fc Hachisul (fTool . we have followed the evolution 
of novae by connecting steady state solutions along the 
decreasing envelope mass sequence. 

The equations of motion, radiative diffusion, and con- 
servation of energy, are solved from the bottom of the 
hydrogen-rich envelope through the photosphere, under 
the condition that the solution goes through a critical 
point of a steady-state wind. The winds are accelerated 
deep inside the photosphere so that they are called "opti- 
call y thick winds." We have used updated OPAL opaci- 
ties (|Iglesias fc Rogers"1ll996[ ) . Optically thick winds stop 
after a large part of the envelope is blown in the wind. As 
one of the boundary conditions for our numerical code, 
we assume that photons are emitted at the photosphere 
as a blackbody with a photospheric temperature Tph. 
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Fig. 6.— Spectra of V838 Her and V1974 Cyg around 1900 
A. In the case of V838 Her, a strong emission line from the C III] 
1907-1909 A doublet is present and appears eroded by overlying 
absorption from the principal and diffuse system components of the 
Fe III UV 34 multiplet. On day 15, these components are seen at 
—550 km and —1600 km s~^, respectively. In V1974 Cyg, the 
C III] doublet emission is not present yet (it will appear at a later 
stage), but it is easy to appreciate the presence of the principal 
and diffuse absorption components of the Fe III UV 34 multiplet 
at -1760 km s'^ and -2930 km s'^ 

The UV 1455 A flux is estimated directly from black- 
body emission instead of modeling WD atmosphere. 
This approximation is reasonably accurate in a very 
narrow wavelength ran ge around UV 1455 A flux 
(jCassatella et al. 1 120021 and paper III) near its peak. 
Infrared and optical fluxes are calculated from free-free 
emission by using physical values of our wind solutions. 
Physical properties o f these wind solutions are given in 
our previous papers (iHachisu et a ."1 Il999al : iKato I Il983l . 
fl997l [l99l [Kato fc Hachisu II1994D . 

We neglect helium ash layer which may develop under- 
neath the hydrogen burning zone in mass increasing WDs 
as m RS Oph (see HachisulFS [2007) . We suppose 
that the WD mass is decreasing in V838 Her, because 
metal enhancement is observed in V838 Her (see Table 
[2]). Metal enrichment is an indicator of decre asing WD 
mass ('e.g.. lPrialnik1ll986l:|Pri"alnik fc Kovetz Ill995[) . 

3.2. Multiwavelength Light Curves 

In the optically thick wind model, a large part of 
the envelope is ejected continuously for a relatively long 
period (e.g.. iKato fc Hachisu I llQM l. After maximum 
expansion of the photosphere, the photospheric radius 
(i?ph) gradually decreases keeping the total luminosity 
(Lph) almost constant. The photospheric temperature 
(Tph) increases with time because Lph = 47ri?pjj(Trp[j. 
The maximum emission shifts from optical to supersoft 



X-ray through UV and extreme ultraviolet. This causes 
the optical luminosity to decrease and the UV luminos- 
ity to increase with a timescale that depends strongly on 
the WD mass and weakly on th e chemical composition 
of the envelope (e.g. iKato Ill997l Paper I). 

3.3. Universal Decline Law of Free-Free Light Curves 

We calculated the optical and IR fluxes of free-free 
emission outside the photosphere of optically thick winds 

by 



K OC / NeN.dV oc 



-f^uh wind 



r dr (x 



-Rph 



(1) 



(see Paper I for more details), where Fi, is the flux at 
the frequency Ne and Ni are the number densities of 
electrons and ions, respectively, V is the emitting vol- 
ume, i?ph is the photospheric radius, Mwind is the wind 
mass loss rate, and Wph is the velocity at the photosphere. 
We also assume that oc Pwind, Ni oc Pwind, and use 
the continuity equation, i.e., pwind — -M^wind/47rr^Wwind, 
where pwind and Wwind are the density and velocity of 
the wind, respectively. Finally, we assume that u^ind = 
const. = Wph outside the photosphere. 

This equation has successfully reproduced the optical 
and infrared light curves of several novae in outburst 
(see papers of t his series; IHachisu fc K"ato1[2rM [200l 
IKato et al.! '20G8). In the case of V838 Her, the early de- 
cay (until day 6) of J, iJ, and K-h&nAs are almost par- 
allel to each other as shown later in Figure [3 whereas 
J is not shown for clarity; This is one of the charac- 
teristic properties of free-free emission (see Appendix of 
Paper HI). Also the spectrum of V838 Her is consistent 
with that of free-free emission as mentioned in Section [TJ 
Thus, we can safely apply our method to V838 Her. 

After the optically thick wind stops, the total mass of 
the ejecta remains constant with time. The flux from 
such homologously expanding ejecta is 



F^<x N^N.dV oc p^V oc -^V oc R 



oc t 



(2) 



where p is the mean density, and Mcj is the ejecta mass. 
We assume that the ejecta are expanding at a constant 
velocity, v. So we have R = vt, where t is the time after 
the outburst. The proportionality constants in equations 
([T]) and ([2]) cannot be determined a priori because radia- 
tive transfer is not calculated outside the photosphere; 
these were determined using the procedure described be- 
low. Fx is also represented by equations ([T]) and (l2|), 
but the proportionality constant depends on wavelength 
since Fx oc \~'^Fy. 

If we assume an expanding shell with constant thick- 
ness, the flux is proportional to F ^ oc as discussed 
in V723 Cas (lEvans et al. l[2003h . In V1500 Cyg and 
V1974 Cyg, however, we found that the flux decline rate 
can be well represented by t~^ (see Paper I). As we will 
see later, the decline rate of V838 Her is consistent with 
t-\ 

As the proportional constant in equation ^ cannot be 
determined theoretically, and there is no representative 
timescale, this later part of the light curve does not con- 
tain any information on the WD mass. Therefore, we do 
not use this part in our light curve fitting. 
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Fig. 7. — UV 1455 A light curve-fitting for V838 Her. Large open 
circles with error bar: observed UV 1455 A flux and root mean 
square of errors (units of the left- hand-side axis). Three different 
sets of chemical composition are assumed, (a) X = 0.55, Xcno = 
0.03, JsTNe = 0.07, and Z = 0.02, (b) X = 0.65, Xcno = 0.03, 
X-Me = 0.03, and Z = 0.02 and (c)X = 0.35, Xcno = 0.2, 
Xf^g = 0.1, and Z = 0.02. The WD mass is attached to each 
light curve. Theoretical flux -F1455 is presented for an arbitrarily 
assumed distance of 1.0 kpc and no absorption (the scale in the 
right-hand-side is for 1.37A/q in (a) and (b), but for 1.35AfQ in 
(c). For other curves scales are normalized arbitrarily). 

3.4. System Parameters of Optically Thick Wind Model 

The light curves of our opticaUy thick wind model are 
parameterized by the WD mass (Mwd), chemical com- 
position of the envelope {X,Y, Xcno, -'^Ne, Z), and the 
envelope mass (AMenv.o) at the time of outburst (JD 
2,448,340.0). Note that the metal abundance, Z, includes 
carbon, nitrogen, oxygen, and neon with solar composi- 
tion ratios and that Xcno and X^c denote additional 
excesses. 

Table [2] summarizes the chemical abundance determi- 
na tions so far available for V83 8 Her, as first provided 
bv lVan landingham et al. ' (1996') , superseded by the later 
re-determinations by Vanlandingham et al. (1997), and 
extended by ISchwarz et al. I (|2007f l from spectra taken 
at four epochs (day 28, 60, 79 and 148). We thus assume 
the chemical composition of V838 Her to be X = 0.55, 
Y = 0.33, XcNO = 0.03, Xno = 0.07, and Z = 0.02. We 
call it the "standard" composition set. We also assumed 
two other sets of composition for comparison. The first 
case is X = 0.65, Y = 0.27, Xqno = 0.03, Xnc = 0.03, 
and Z — 0.02, expressing less enrichment of heavy ele- 
ments such as in V382 Vel and QU Vul, and the second 
case ia X = 0.35, Y = 0.33, Xcno = 0.2, = 0.1, 
and Z — 0.02 which may be an extreme case of metal 
enrichment for V838 Her but appropriate composition 
for V1974 Cyg and V351 Pup, although abundance es- 



timates are very different among different authors (see 
Table 1 in Paper I which lists observational estimates of 
abundance of nova ejecta). 

The latter two recent determinations both indicate 
an excess of h elium with respect to hydrogen, i.e., 
X/Y = 1.9 ([Vanlandingham et aH [T997il and 1.8 
(jSchwarz et al. II2007D . against the solar value of 0.7/0.28 
= 2.58. In massive WDs about one tenth of hydro- 
gen is consumed by nuclear burning to produce ther- 
mal and gravi tational energy during the early phase of 
the o utburst (|Politano et al. I [l995 : Pri alnik fc Kovetzl 
Il995f ). Supposing that AX = 0.1 of hydrogen is con- 
verted into helium, the hydrogen/helium ratio becomes 
X/Y = (0.7 - 0.1)7(0.28 + 0.1) = 1.6. If, in addi- 
tion, a fraction 10% of the envelope mass is dredged up 
and mixed into the envelope, the above envelope com- 
position is further changed to X ^ 0.545, Y = 0.345, 
and XcNO + -'^Ne + Z = 0.11. This is very consistent 
with t he o bservational resu l ts by IVanlandingham et al. I 
(IT991 and ISchwarz et al. I (l2007l ) in Table H 

4. LIGHT CURVE FITTING 

We have searched for the model light curves that best 
fit the observational data among models with WD mass 
1.2, 1.3, 1.33, 1.35, and 1.37 Mq, for the three sets of 
chemical abundance quoted above. As shown in previ- 
ous papers of this series, we have proven the light curve 
fitting method to be a powerful tool to obtain the intrin- 
sic nova parameters. In particular, we have shown that 
the theoretical light curves of UV 1455 A continuum are 
very sensitive to the WD mass, so that the WD mass can 
be determined with a good accuracy. 

Figure [7] shows the fits to the UV 1455 A hght curve of 
V838 Her for three different chemical composition and 
various WD masses. Figure [7^) shows that Mwd = 
1.37 Mq model is in a good agreement with the observed 
data but models of Mwd < 1-33 Mq produce too slow 
declines. Models of 1.2 and 1.3 Mq are omitted because 
they are too slow to be compatible with the observa- 
tion. Similar results are obtained if we increase the hy- 
drogen content and decrease the heavy elements slightly 
to (X,XcNO,>^Ne) =(0.65,0.03,0.03) as in Figure [7)3. 
On the other hand, if we decrease the hydrogen and in- 
crease the heavy elements to X = 0.35, Xcno = 0.2, and 
-'^Ne = 0.1 (Figure[7t), which may be an extreme case for 
V838 Her, the 1.35M0 model shows a good fitting and 
less massive WDs (Mwd < 1-3) are probably excluded. 
Therefore, from the UV light curve fitting, we see that 
the WD in V838 Her is as massive as Mwd > 1-35 Mq. 

Figure [8] reports the model fits to the visual and in- 
frared light curves for the standard composition set. We 
have assumed that, at these wavelengths, the light curve 
is dominated by free- free emission (i.e., can be calculated 
from eq. [1]). Since its shape does not depend on wave- 
length, one can apply the same light curve model to the 
optical and IR data. Figure [5] shows indeed that, if suit- 
ably upward-shifted, the optical light curve also fits the 
infrared fluxes, at least until day 6, when an IR flux min- 
imum was reached (indicated as a vertical dashed line in 
Figure [5]), owing to thermal emission from dust (see Sec- 
tion [T] and Section [5|). 

Figure [9] provides, in logarithmic time, another view of 
the model fit to the data. This figure shows that, at very 
late phases, optical and infrared fluxes decay along with 
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Fig. 8. — Theoretical and observed optical, infrared, and UV 
light curves of V838 Her. Theoretical light curves refer to models 
with AfwD = 1-37 Mq (blue), 1.35 Mq (black), and 1.33 Mq 
(red). Observed visual and V magnitudes are denoted by small 
open circles (lAU Circulars 5222, 5223, 5224, 5226, 5228, 5229, 
5233, 5234, 5238 , 5254 and 5265), and by diamonds {lUE VpES- 
ICassatella et al. I l l2004al ')'). Two large downward arrows indicate 
optical upper limit observations from lAU Circular 5265. Small 
downwa rd arrow s denot e the da tes of eclipse mi n ima r eported by 
ILeibowitz eFain | |T993) and by |Kato & Hirata I flM^. Infrared 
data are taken from lAU Ci rcular s 5223, 5224, 5228, 5230, 5232, 
5236, 5246, 5254 and from [Harrison fc Stringfellowl (HQQI 'I. K 
and H magnitudes are denoted by middle-sized open circles and 
crosses, respectively. H magnitudes are shifted up by 0.5 mag. 
Large open circles denote the observed UV 1455 A band flux 
iCassatella et al. ,2002 . except one newly obtained point) in units 
of 10"^-^ erg s~^cm~'^ A~^. For theoretical UV curves, we assume 
the distance and absorption in Table|3]to fit the observational data. 
Vertical dotted line indicates the epoch of the IR minimum. 

the Fa oc line as predicted from equation ([2]). This 
circumstance does not bring, however, any information 
on the WD mass. As shown in Figures [H] and [51 among 
models with the WD mass of 1.33 , 1.35 and 1.37 Mq, 
the 1.35 Mq WD provides a best fit representation of 
the UV light curves as well as the optical and infrared 
light curves. 

Figures [TUl and [TT] show the model fits for the other two 
sets of chemical composition. In Figure [TOl we get the 
best fit models of 1.35 Mq WD from fitting optical, IR 
and UV fluxes. In FigurefTT] however, we cannot fit both 
of the UV and optical fluxes simultaneously. When we fit 
the optical, the model UV flux rises too earlier, and then 
we cannot obtain a good fit in the UV. This indicates that 
the assumed composition, i.e, X — 0.35, Xcno — 0.2, 
and Xnc = 0.1 is too metal rich and not appropriate for 
V838 Her. 

To summarize we may conclude that the WD mass is 
M-yvD = 1-35 ± 0.02 Mq. The model parameters and the 
main results are summarized in Table [31 

5. DISCUSSION 

5.1. "Optical Drop" in \ Woodward et al. I U99^ ) 

[Woodward et alH ()1992[ ) and [Harrison fc Stringfellowl 
(Il99l " claimed that the rise of infrared fluxes starting on 
JD 2,448,346 (day 6) was accompanie d by a sudden drop 
of the visual magnitude by AV^ ~ 0.9. [Woodward et al. 1 
(|1992l ) proposed that the sudden drop of the visual flux 
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Fig. 9. — Same as Figure [S] but for logarithmic time. The 
AAVSO data (small dots) are plotted only until day 68, be- 
cause after that time the nova becomes fainter than 14 mag, i.e., 
the brightness of a closeby star. The optical data f or the late 
phases (filled squares) are taken from llngram et al. I l|1992l ) and 
ISzkodv fc In gram (1994). The IR data denoted by the asterisk 
are taken from Harrison fc StringfcUow (1994). The emergence 
time of companion for the 1.35 Mq WD model (day 9) and the 
first eclipse observation (day 21) are also indicated by the small 
and large downward arrows. 
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Fig. 10. — Same as Figure |9] but for theoretical model with 
the composition set of less enhancement of heavy elements, i.e., 
X = 0.65, XcNO = 0.03, and = 0.03. 

was caused by an optically thick clump of matter ejected 
toward the earth on day 6 (indicated by the vertical dot- 
ted line in Figure [12]). Afterwards, the optical flux con- 
tinued to decline but at a substantially lower level than 
extrapolated from the light curve of previous days, shown 
by a dash-dotted line in F igure | 12 | (which mimics the 
dashed curve in Figure 1 of [Woodrord et al. l[T99l. 



The visual m agnitudes used bv [Woodward et al. I 



(|1992f) and by [Harrison fc Stringfellow [ (|1994[ ) were 
taken from lAU Circulars. These data were, however, 
obtained by different observers. We have tested self- 
consistency of the different sets by connecting the data 
obtained by the same observer (Figure [T^ . It clearly 
appears from the figure that there are systematic dif- 
ferences between the data sets from different observers, 
and also that no observer reported the sudden drop at JD 
2,448,346, in agreement with the smooth decline shown 
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TABLE 3 

Summary of the present model 



subject units model 1 model 2 model 3 model 4 model 5 

1.37 Mq 1.35 Mq 1.33 Mq 1.35 Mq 1.33 Mq 



outburst day (to) 




JD 2,448,340 


E{B - V) 




0.53 


X 




0.55 


Y 




0.33 


CNO 




0.03 


Ne 




0.07 


Z 




0.02 


distance from UV fit 


kpc 


2.7 


My^niax from UV fit 


mag 


-8.41 


WD envelope mass 


■ Mq 


2.5 X 10-^ 


mass lost by winds 


■ Mq 


2.3 X lO"*^ 


wind phase 


days 


37 


secondary mass 


■ Mq 


0.76'' 


separation 


■ Rq 


2.4 


Rl 


Rq 


1.0 


Ri 


. Rq 


0.8 


companion's emergence 


days 


10 



^ estimated from cq. JlOt . 




12 3 



log[days after outburst] 

Fig. 11. — Same as Figure |9] but for theoretical model with 
the composition set of much enhancement of heavy elements, i.e., 
X = 0.35, XcNO = 0.2, and Xn,, = 0.1. 

by the homogeneous set of the lUE Vfes magnitudes, 
also given in the same figure. Moreover, the magnitudes 
corresponding to just before and after the sudden drop 
are 8.5 mag at Mar. 30.46 UT and 9.54 mag at Mar. 
30.466 UT, respectively, which are obtained by different 
observers. If we take the reported observation time as 
a three-digit number, the magnitude should dropped by 
0.006 day — 8.6 minutes. This is too short to have some 
significant change in stellar luminosity. Therefore, we 
conclude that the claimed V magnitude drop is an ar- 
tifact caused by connecting the different data set with 
systematic difference. 

Should the 0.9 mag optical drop be real and due to 
dust absorption from a clump of matter ejected along the 
line of sight, one would expect to observe a simultaneous 
strong absorption of the UV flux. The decrease in the UV 
flux is estimated to be a factor of lo(8-3/3 i)*o.9 = 257. 
Figure [8] shows, however, no indication of such a drastic 
decrease in the UV flux. 

5.2. Dust Formation: Comparison with V1668 Cyg 

After the local minimum on day 6, the IR fluxes 
increased to reach a local maximum around day 12 
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Fig. 12. — Optical and infrared light curves. The optical data are 
shown separately to assess their homogeneity. From bottom to top, 
IAUC-l-3: data from lAU Circulars, shifted by 3 mag downward, 
in which all the data are connected along the reported observing 
time. lAUC: same data of lAU Circulars but are connected be- 
tween points only obtained by the same observer. lUE Vfes — 2: 
lUE VpES data connected by a dashed line (shifted 2 mag up- 
ward). K — 1 and H — 1.5: infrared K and H data, shifted by 1.0 
and 1.5 mag upward, respectivel y. The time of th e "drop in optical 
magnitude" proposed by Woodw ard et al~l l|1992l ') is indicated with 
a vertical dashed line, and their "unabsorbed light curve" with a 
dot-dashed line. 

(in H band), being ascribe d to thermal emission from 
an optically thin dust shell fChandrasekhar et al.''1993; 
IWoodward et al. 1992; Kidger & Martinez-Roger 1993 
[Harrison fc Stringfellow 1 11994[ ) . As shown in Figure [H 
however, there is no clear evidence for a corresponding 
decrease in the UV fluxes. In this subsection we make a 
comparison with nova V1668 Cyg in order to understand 
to what extent dust absorption can cause significant ef- 
fects on the UV and optical fiuxes in the case of V838 
Her. 

V1668 Cyg is a well studied moderately fast nova. It 
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developed an optically thin dust shell giving rise to the 
infrar ed second maxirnum ^ 60 days after the optical 
peak (jStickland et al. Ill981h . In correspondence to the 
second IR maximum, the UV 1455 A flux suddenly 
dropped by 1/2 ^ 1/3 of what exp ected from a linear 
decHne law (jKato fc Hachisu Il2007f l. The evolution of 
V1668 Cyg (^3 ~ 23 days) was much slower than that 
of V838 Her, due to its smaller WD mass 0.95 Mq, 
see Paper I). Consequently, the ejected mass of V1668 
Cyg is expected to be much larger than that of V838 
Her. Namely, since the ejecte d mass of V1668 Cyg i s 
estimated to be 5.8 x 10"^ Mq (jKato fc Hachisu Il2007t ). 
while we estimate that it was 2 - 3 x 10"^ M© for V838 
Her from our best fit models (see Table [3]). One should 
thus expect that the dust mass of V838 Her was roughly 
ten times smaller. 

Observationally, iGehrz et al. I (|1980( ) estimated for 
V1668 Cyg the mass of dust (condensed carbon) to 
be 2.5 X 10~^ Mq and shell mass of hydrogen to be 
2 X 10~^ Mq, which is consistent with the above value 
of total ejecta mass 5.8 x 10~^ Mq 

Also, IChandrasekhar et al.l (|1992( ) estimated the total 
grain mass of V838 Her to be 1.6 x 10~^ Mq assuming 
distance of d — 8.3 kpc. As the mass estimate is propor- 
tional to d^, where d is the distance to the star, this cor- 
responds to 1.7 X 10~^ Mq for d = 2.7 kpc here adopted 
(see Section 5.3 below). iHarrison fc StringfelloVI (fT99i ) 
estimated the dust mass of V838 Her to be 1.1 x 10~* Mq 
at dust shell maximum for d — 5.5 kpc, which is propor- 
tional to the IR luminosity, therefore, a smaller value 
2.7 X 10-^ Mq is derived for d = 2.7 kpc. 

In any case, the dust mass of V838 Her appears to be 
one tenth of dust mass of V1668 Cyg, a value that is not 
expected to affect significantly the UV fluxes even at the 
time of the second infrared maximum, explaining then 
the absence of any quick drop in the UV 1455 A light 
curve. 

5.3. Distance 

The interstellar extinction and the distance were esti- 
mated with__vaji£u^^ 

119941: IVanlandingham et al. 1 1 19961 for summary), the 
values of which are largely scattered as E{B — V) = 
0.3 - 0.6 and d ^ 2.8 ~ 10 kpc (see Table HJ. One of the 
typical way to estimate the distance is to use ^3 time or 
using statistical properties such as a given magnitude at 
a certain stage of light curve. 

Let us begin with the distance estimate using the Max- 
imum Magnitude versus Rate of Decline (MMRD) law. 
The absolute V magnitude at maximum Afy.max can 
be estimated fror n the ^3 time th rough Schmidt-Kaler's 
MMRD relation (jSchmidt \\195% : 

Mv,max = -11.75 + 2.5 log i3. (3) 
Then the distance is obtained from 
(m - M)v,max = -5 + 51og d + 3.1EiB-V), (4) 

if E{B — V) is known. If we adopt ^3 = 5.3 from the 
optical maximum of my = 5.4 ()Harrison fc St ringfcllow ' 
[1994) . this equation gives My = —9.94. With the appar- 
ent visual magnitude of my — 5.4 and E{B — V) — 0.53, 
we get d = 5.49 kpc. If E{B — V) is not fixed, we have 

(TO-M)vmax = -5 5 log d 3.1E{B-V) = 15.34, 

(5) 




0.5 1 



E(B-V) 

Fig. 13. — Distance-reddening relations for V838 Her. Solid 
line labeled "UV 1455 A" (see Eq. [8]), corresponds to our best 
fit model MwD = 1-35 Mq. Dotted curve represents the max- 
imum magnitude vs. rate of decline relation calculated from 
Schmidt-Kaler's law, with ts = 5.3 days (labeled "MMRDl", eq. 
\5\ ) ; the dash-dotted cu rve shows another MMRD relation, by 
IDella Valle fc Livio I ilMB ). with (2 = 1.4 day (labeled "MMRD2", 
eq. [7]). "BB" is the line of eq. l[9]l. An arrow indicates our redden- 
ing determination E{B — V) = 0.53. Open circles denote ne arby 
stars' distance vs. E{B — V) (taken from (46/2) iniNeckel & Kl are I 
iflQSOl) ). See text for more details. 

which is plotted in Figure [H (labeled "MMRDl"). 

The Schmidt -Kaler law ( Eq. ([51)) yields a bit larger 
magnitudes for a faster nova. IDella Valle fc Livio I (fl995h 
derived the following MMRD relation which is a good 
representation of fast novae, 

Mv,max = -7.92 - O.Slarctan ^'^^"^"^^^'^ (6) 

For t2 = 1.4 dH arrison fc Stringfello"w1fT99l this equa- 
tion gives My = —9.04. With the apparent visual mag- 
nitude of my — 5.4 and E{B — V) = 0.53, we get d = 3.6 
kpc. If E(B — V) is not fixed, we have 

(m - M)v „iax - -5 + 5 log rf + 3.1E{B -V) = U.U, 

(7) 

which is plotted in Figure [13] (labeled "MMRD2"). 

Another way to estimate the distance to V838 Her 
comes from comparing the observed light curve of 
1455 A band with the corr esponding model fluxes 
(see Paper I, Paper HI, and iKato fc Hachisu I 120051 
I2007f ). We show that the calculated flux at A = 1455 
A at a distance of 1 kpc for Model 2 is F^°'^ = 6.6 
xlO"-^" erg cm~^ s~^ A~-^ at the UV peak in Figure 
[7^. The corresponding observed peak flux is F^'^*'=1.57 
xlO^^^ erg cm~^ s~^ A"-*^. From these values one ob- 
tains the following relation between distance and redden- 
ing: 

^obs _ ^^mod ^ 5 log(^) + AxEiB -V)^ 6.56 (8) 
1 kpc 

whe re mx = -2.5 1og(FA), and Ax = 8.3 at A = 1455 
A ([Seaton I I1979D . This curve is reported in Figure 
[H and is labeled "UV 1455A". Equation ^ provides 
a distance of 2.7 ± 0.5 kpc for E{B - V) ^ 0.53 ± 0.05. 
If we take the upper and lower values of the UV peak 
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(1. 57±0. 20) X 10^12 ergcm-^g-i we get the dis- 

tance of 2.9 kpc and 2.5 kpc, respectively. This is how- 
ever, an extreme case, because the fitting becomes much 
worse in Figure [7l Therefore, we estimate the distance 
to be 2.7 ± 0.5 kpc for E{B ~V)^ 0.53 ± 0.05 from UV 
light curve fitting. Note that the distance hardly changes 
for all the models in Table [31 

One may estimate a lower limit for the distance assum- 
ing that the flux of free-free emission is larger than the 
blackbody flux in optical wavelength region. As the op- 
tical free-free flux is dominated by photons reprocessed 
from shorter wavelength photons of blackbody contin- 
uum, free-free flux is larger than that of corresponding 
blackbody emission at the opti cal wavelength during th e 
early phase of nova outburst fHauschil dt et ainil997[ ). 
Our experience of light curve fitting also suggest that 
this condition is empiricall y satisfied in several novae 
Csee iKato fc Hachisu I [2007L for V1974 Cyg). Thus, we 
assume this condition, which is written as 

TOv -Mv = 51og(— ^) + 3.1E(B -V) > 13.2. (9) 
10 pc 

The shadowed region labeled "BB" in Figure [Ol is where 
Equation ^ is not satisfied. The border gives the lower 
limit for the distance to V838 Her, which is d > 2.46 kpc 
for the particular value of E{B — V) = 0.53. 

As shown in Table (TJ the distances in literature are 
very scattered. The distance obtained from MMRD 
tends to be larger th an ours as shown in Figure 1131 
iStarrfield et al. I (| 1992) obtained a distance of ^ 3.4 kpc 
and E{B — V) ^ 0.6 from the comparison of total UV 
flux (LWP/SPW) including both of line and continuum 
emission with that of nova LMC 1991 showing a similar 
spectrum development. This is similar but contrasted 
to our method in the sense that we use only continuum 
emission around 1455 A and compare it with the theoret- 
ical value. Their value ^ 3.4 kpc is larger than ours, but 
closer than other estimates using the MMRD relation. 

5.4. Emergence of the Companion 

In the early phase of the outburst, the photosphere of 
the WD envelope extends beyond the size of the binary 
orbit. The companion is deeply embedded inside the 
WD photosphere. After maximum expansion, the photo- 
spheric radius shrinks owing to wind mass loss. We can 
estimate the epoch when the companion appears from 
the photosphere. 

From Warner's (1995) empirical formula the mass of 
the donor star (companion) is expressed as: 

^ « 0.065 ( , for 1.3 < < 9 (10) 

Mq \ hours / hours 

which gives Ma = 0.76 M© for Forb = 0.2976 days 
(7.14 hr), in good agreeme nt with the value M2 = 
0.73 - 0.75 Mq obtained bv iSzko dv fc Ingram I ([1991 . 
although they could not well constrain the mass of the 
primary component (AfwD > 0.62 Mq). 

The binary separation is obtained from Kepler's law 
to be a = 2.41 i?© for Mwd = 1-35 Mq. The effective 
radius of the Roche lobe is Rl — 1.0 Rq for the primary 
component (WD), and i?2 — 0.80 Rq for the secondary. 
Here we u s e an e mpirical formula of Rl and R2 given by 
lEggletonl (|1983D . In our model, the companion emerges 



from the WD envelope when the photospheric radius of 
the WD shrinks to Rph ~ 3.2 Rq (the separation plus 
R2), ~ 2.4 Rq (the separation), or 1.6 Rq (the sepa- 
ration minus R2). This occurs on day 7, 9 and 11, re- 
spectively, in our Model 2. The emergence time of the 
companion (the central time) is shown in Table [3l and 
also indicated by a small arrow in Figures M [IHI and [TlJ 
As the eclipse is due to the occultation of the WD, or of 
the central part of the accretion disk by the companion, 
the eclipse must occur after the above three epochs of 
emergence. Our estimated time of companion emergence 
is a rough estimate, but it is well before the first eclipse 
observation (day 21), so consistent with observation. 

5.5. Mass determination from the duration of the 1455 
A maximum 

The duration of a UV nova outburst, as measured 
by the full wid th at half maximum ^f whm of the 1455 
A light curve (jCassatella et al. |[20M ) is tightly corre- 
lated with the nova decay time ^3 which, in turn, is a pri- 
mary indicator of the WD mass (Livio 1992). Detailed 
model calculations bv lHachisu fc Kato I (|2006D have con- 
firmed the primary dependence of tpwHM on Mwd and 
have shown, in addition, that there is secondary, weaker, 
dependence on chemical composition. These results im- 
ply that the duration of a UV nova outburst can be used 
as a useful tool to estimate the WD mass in nova systems. 

In Figure [TH we comp are the theoretica l relation of 
^FWHM vs. WD mass (jHachisu fc Kato I [2006) with 
those obt ained from t he measured values of tpwHM 
(Cassatcll a et al. I l2002l ) and from the WD masses of 
seven classical novae, including V838 Her, estimated 
from the previous works (na mely from Paper I for V I 668 
Cyg and V1974 Cyg; from IKato fc Hachis^ [200I for 
V351 Pup, OS And and V693 CrA, and from this paper 
for V838 Her) . We recall that the masses of these novae 
were determined by using the light-curve fitting method 
not only on the UV 1455 A light curve, but also on the 
optical, infrared, and supersoft X-ray light curves. For 
GQ Mus the UV data are very scattered so we plot here 
the theoretically fitted value in Paper HI. 

The largest possible error in this mass determination 
arises from ambiguity of chemical composition because 
observationally determined values are sometimes very 
scattered in literature. If we adopt different chemical 
composition, we get different AfyvD even for the same 
^FWHM as shown in Figure [T4l In the case of V838 Her, 
however, the difference of mass determination AMwd is 
relatively small, because in very massive WDs, ^fwrm 
quickly decreases and then strongly depends on the WD 
mass but weakly on the composition. 

Another possible source of ambiguity in the determina- 
tion of WD mass may come from effects of dust. Figure 
[T5|compares the duration of UV nova outburst tpwHM, as 
predicted by our models, with the observed values from 
Cassatella et al. (2002) for five novae. The figure shows 
that there is a reasonable agreement between the two sets 
of data, except the dust forming novae V1668 Cyg and 
OS And. In these two objects, dust formation took place 
shortly after the UV 1455 A maximum, so making the 
measured ipwHM significantly smaller than expected. In 
the case of V838 Her, in spite of being a dust forming 
nova, the mass of dust was too small to affect the UV 
light curve, as discussed in Section [5?2l 
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Fig. 14.— Duration (FWHM) of UV nova outburst at 1455 
A against the WD mass. Dash-dotted line: Neon novae {X = 
0.65, Y = 0.27, XcNO = 0.03, = 0.03). Top solid line: Neon 

novae {X = 0.55, Y = 0.30, Xcno = 0.10, Xno = 0.03). Dashed 
line: CO novae {X = 0.55, Y = 0.23, Xcno = 0.20). Bottom solid 
line: CO novae (X = 0.45, Y = 0.18, Xcno = 0.35). Dotted line: 
CO novae {X = 0.35, Y = 0.33, Xcno = 0.30) (almost overlapped 
with the lower solid line). We assume Z = 0.02 for all the models. 
Arrows indicate the WE) mass of individual novae. For V838 Her 
we adopt Model 2. For GQ Mus, we used the results in Paper 
III. For V1668 Cyg and V1974 Cyg, the result s in Paper I. For 
¥351 Pup, OS And and V693 CrA, the results in lKato fc Hachisu I 
1120071 ). Novae named in the upper side of the lines are Ne novae 
whereas those in the lower side are CO novae. 
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Fig. 15. — Observational and theoretical durations of the UV 
1455 A nova outbursts. Filled symbols (red) denote novae with 
dust formation. Squares and circles denote Ne and C O novae, 
respec tively. Observational data are taken from Cassatell a et al.~l 
mo^ . The large deviation of OS And and V1668 Cyg is due to 
dust formation. See text for details. 



Note that the mass determination of our light curve 
analysis is based not only on the UV data fittings, but 
also on the optical and IR fittings, so the error in UV 
fitting docs not directly refiect to the resultant Mwd- 

Figure [Til indicates that less massive WDs (< 1 Mq) 
are systematically CO-rich, while more massive WDs (> 
1 Mq) are neon-rich. It has been argued that heavy 



element enrichment of nova ejecta is caused by dredge- 
up of WD core material so that, for example, C and O 
enrichment is an indication of a CO WD. An ONe WD 
is born in massive stars after ignition of a CO core and, 
therefore, ONe WDs are usually more massive than CO 
WDs when they are just born. Our results are consistent 
with this. 

The mass boundary of CO WDs and O Ne WDs are 
suggested to be 1.0-1.1 Mq. Ilben I (|199ClD summarized 
that a star of initial mass of 5 -10 Mq leaves a CO 
WD of mass 0.75 - I.IMq and a star of 10 - 12Mq be- 
comes ONe WD. iGarci'a-Berro et al. I (|1997[ ) calculated 
stellar evolution of 9 Mq star which finally becomes an 
^ I.OSMq ONe core surrounded by an CO envelope of 
^ O.OSMq. lUmeda et aTl (|1999) calculated stellar evo- 
lution until CO core ignition with fine grid of initial stel- 
lar mass and concluded that the mass of CO WDs born 
in binaries is less massive than 1.0 7 Mq of wh i ch ini - 
tial mass is - 8Mq for Z = 0.02. iMeng et al.~l ()2008f ) 
obtained CO core mass consistent with Umeda et al.'s 
value. These predictions are consistent with our results 
in Figure [Ml 

5.6. X-ray observations 

X-ray emission from V838 Her has been detected 
with ROSAT as early as 5 days afte r the discovery 
(|Llovd et al. l[l99l[OWn et al. I[l99l . and was found 
to be consiste nt with emission from shocked plasma with 
r 10 keV (jO'brien et al. Ill994f) . According to these 
authors the shock could originate from within the ejected 
material itself or by interaction of the ejecta with pre- 
existing circumstellar matter. In our Model 2 the opti- 
cally thick wind lasts until day ~ 52 days (see Table [3] 
for other models); this is consistent with the picture that 
the wind collides with circumstellar matter to produce 
the observed X-ray fiux. 

5.7. Quiescent Phase 

The presence of an accretion disk has been suggested 
from the long duration of the eclipse (jlngram et al. I 
1 1 9921). and from the depth of eclipse at minimum 



( Leibowitzlll993D . 

The absolute magnitude of a disk, seen at an inclina- 
tion angle is approximated by 



Mv(obs) = -9. 



5 / Mwd Mace 

3 Mq Mq yr-1^ 



--log(2cosi), 



(11) 



IS the WD mass, and M^^.r. is the mass 



where M-wd 

accretion rate (E g. [A6) in IWebbink et al.l Il987f) . 

ISzkodv &: Ingram"! (|1994D estimated the inclination of 
the disk to be i = 78 - 90°. Assuming Mwd = 1-35 Mq 
and i = 80°, we get My = 6.5, 4.8, and 3.1 for 
the mass accretion rates of 1 x 10^^, 1 x 10~*, and 
1 X 10~^ Mq yr~^, respectively. The corresponding ap- 
parent magnitudes, as calculated from equation (j4|) with 
E{B -V) = 0.53 and d = 2.7 kpc are my = 20.3, 18.6 
and 16.9. Considering the ambiguity of assumed accre- 
tion rates and the inclination angle, these values are con- 
sistent with the preoutburst magnitude reported in Sec- 
tion [H 
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The contribution of a companion st ar to the qu i escent 
luminosity is estimated as follows. iPattersonl (|1984l ) 
gives an empirical relation between the absolute mag- 
nitude of a Roche lobe filling companion and its orbital 
period as 

My ^ 22 - 17.46 logP(hr), (12) 

for 0.7 < logP(hr) < 1.1 . For -P(hr) = 7.14 
hr ([Ingram et al. 1119921 : iLeibowitz et al. I I1992D we get 
My = 7.1 and my = 20.6. Thus, the companion star 
would be fainter than the accretion disk although there 
are ambiguity in the mass accretion rate and the incli- 
nation angle. This is also consistent with a faint com- 
panion suggested from a small upper limit of ^ 0.05 
mag on the depth of the secondary minimum (|Leibowitz I 
Il993f) and from no visi ble evidence in the spectrum 
(jSzkodv fc Ingram |[l99l . 

6. CONCLUSIONS 

We have applied the 'universal decline law' of classical 
novae described in §3.31 to one of the fastest novae V838 
Her and derived various parameters. Our main results 



are summarized as follows: 

1. An analysis of the lUE reprocessed data of V838 Her 
indicates E{B -V) = 0.53 ± 0.05. 

2. The 1.35±0.02 Mq WD model reasonably reproduces 
the light curves of V838 Her both in the optical and in 
the UV 1455 A band as well as in the infrared H and 
K bands at the early stages, until dust formation takes 
place. Model 2 in table [3] is the best fit model. 

3. The distance is estimated to be d ~ 2.7 ±0.5 kpc from 
the UV 1455 A hght curve fitting. 

4. We have estimated ejecta mass AMwind '--^ (2 — 3) x 
1O~'^M0 lost by winds. 

The authors are grateful to Takashi lijima for fruitful 
discussion, and also to an anonymous referee for useful 
comments to improve the manuscript. We also thank 
the American Association of Variable Star Observers 
(AAVSO) for the visual data of V838 Her. This research 
has been supported in part by the Grant-in- Aid for Sci- 
entific Research (20540227) of the Japan Society for the 
Promotion of Science. 



APPENDIX 

Figure [16] demonstrates a "universality" of nova light curves, in which five classical novae are shown in the visual 
{V and y) and UV 1455 A band (except V1500 Cyg: no UV 1455 A band is available). Each peak of the UV 1455 
A flux is normalized in order to match the peak of V838 Her. Also, the time-scale of each nova is normalized so that 
each UV 1455 A light curve is overlapped with that of V838 Her. The same normalizing factor of time-scale is applied 
to both the UV and the optical light curves. These normalizing factors, obtained by eye fitting, are summarized in 
the figure caption. 

Our universal decline law predicts that the light curves, normalized in this way, should merge into one, independent 
of the WD mass and chemical composition (see Pape r I). This match is e xcellent for the y light curve, because this 
band is little contaminated by strong emission lines l[Hachisu et al. ir2008| ). As one may appreciate from the figure, 
there is a good overlap between the light curves of V1668 Cyg (y-magnitude, connected by solid line), V1500 Cyg 
(y-magnitude) , and the initial phase of V1974 Cyg {V magnitude). It also appears from the figure that, whenever 
strong emission lines provide the dominant contribution to the magnitude, the data deviate above the V1668 Cyg line. 
In the case of V838 Her the data are scattered compared with the above three objects, therefore, we fit its lowest 
boundary of the data with the V1668 Cyg line. We can see excess of the visual magnitudes between day 3 and day 6. 

Four UV 1455 A light curves also show a good agreement with ea ch other, although V1668 Cyg shows a drop at 
0.37 in the normalized time due to dust formation ([Gehrz et al. I[l980l ). We don't see any indication of sharp drop like 
this in V838 Her. 
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Fig. 16. — Scaling of five nova Ught curves. The observed Ught curve of classical novae, V693 CrA, V1974 Cyg, V1668 Cyg, and V1500 
Cyg are overlapped to that of V838 Her. For V838 Her only the optical data in lAU Circulars and lUE VpES data are shown. Optical data 
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